Testing Quintessence models with large-scale structure growth 
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We explore the possibility of putting constraints on quintessence models with large-scale structure 
observations. In particular we compute the linear and second order growth rate of the fluctuations 
in different flavors of quintessence scenarios. We show that effective models of quintessence (e.g. 
with a constant equation of state) do not account of the results found in more realistic scenarios. 
The impact of these results on observational quantities such as the shape of the non-linear power 
spectrum in weak lensing surveys or the skewness of the convergence field are investigated. It appears 
that the observational signature of quintessence models are specific and rather large. They clearly 
cannot be mistaken with a change of Qo- 
Pacs numbers: 

98.80.Cq, 98.65.Dx, 98.62.Sb 
I. INTRODUCTION 

The recent evidences in favor of a non-zero cosmologi- 
cal constant Jl|-§| have led to developments of alternative 
scenarios for explaining such a non-zero vacuum energy 
density pPf . In particular the models involving the so- 
called "Quintessence" have attracted attention from the 
high-energy physics community |11-17|. Indeed in such 



an unavoidable parameter degeneracy [|21],[22| so that 
cannot be determined alone. 

The impact of quintessence models on the properties 
of the CMB as well as the primordial density contrast 
power spectrum has nonetheless been studied in different 
models: the effective Quintessence (with the extra short- 
coming that the possible fluctuations of the quintessence 
field were neglected) as well as some high energy physics 



models the vacuum energy density is due to the poten- 
tial and kinetic energy of a scalar field rolling down its 
potential. Various models has been proposed. The sim- 
plest implementation of such models, widely used in the 
literature, is to introduce an effective quintessence with a 
constant equation of state |fl). More elaborate theories 
provide potentials that exhibit a tracking solution regime 
as long as the energy density of the quintessence field 
is sub-dominant |p_3|| . This behavior is generically en- 
countered in the Ratra-Peebles |f7| model in which the 
quintessence potential is a simple inverse power of the 
field. Other models inspired from high-energy physics 
have also been shown to exhibit such a remarkable prop- 
erty 

The presence of a quintessence field changes the en- 
ergy content of the universe and therefore alters its global 
expansion rate. It is then natural to try to detect the 
signature of a non-standard vacuum equation of state 
through its impact upon the distance-luminosit y f unc- 
tion as it can be revealed by SNIa observations Q po| . 
It has been found however that the precision with which 
the vacuum equation of state can be measured depends 
crucially on whether priors are assumed on the other cos- 
mological parameters in particular on the matter content 
of the Universe. It calls for a re-examination of the theo- 
retical foundations upon which precision methods for the 
determination the cosmic density are based. 

The Cosmic Microwave Background anisotropy power 
spectrum has been recognized as a gold mine for the de- 
termination of the cosmological parameters. It is actually 
a very precious way for measuring the global curvature 
of the Universe |lj (through the value of the angular dis- 
tance of the last scattering surface) but it suffers from 



tracking potential [14 231]. It has been found that at the 



redshift of recombination the dark energy fluid is sub- 
dominant and has only significant super-horizon fluctu- 
ations. Quintessence effects reveal therefore only as a 
modest change of the Sachs- Wolfe plateau, an effect dif- 
ficult to be unambiguously detected because of the im- 
portance of the cosmic variance. 

It has been shown however that if the intrinsic proper- 
ties of the CMB anisotropies fail to provide for an unam- 
biguous test of quintessence, its existence can be betrayed 
by the amplitude of the density fluctuations on the last 
scattering surface compared to those at low redshift. This 
can be done for instance with the help of galaxy cluster 
counts |24|,^5| or with weak lensing measurements |^6],^7j . 

In all cases however direct constraints of VLq that would 
help to disentangle models rely on analysis of the local 
universe properties. Matter content of galaxy clusters or 
their number density evolution |24|,|2^] can provide use- 
ful constraints. Unfortunately these methods depend on 
non-trivial modeling of cluster properties such as X-ray 
luminosity or temperature-mass relations. It is therefore 
unlikely they can provide accurate constraints on Qq with 
a well controlled level of systematics. 

New methods, based on weak lensing observations, are 
now emerging that are in principle free of elaborate phys- 
ical modeling. The proposed means for constraining Qq 
are based on the rate at which non-linear effects start 
to play a role in the cosmic density field. Fundamen- 
tally two ideas have been followed. One is based on the 
nonlinear evolution of the shape of the power spectrum 
p9| and preliminary results have already been reported 
in this case ^7|. It relies here on some specific class of 
models namely some flavor of CDM model (the shape 
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of the nonlinear power spectrum depends obviously on 
what is assumed for the linear one). The other method, 
proposed in |30| , is more demanding on the observation 
side but is based on the only assumption that the initial 
conditions were Gaussian. It relies on the direct detec- 
tions of non-Gaussian properties of the density field. In 
particular it has been shown that the large-scale conver- 
gence skewness can be used to measure Hq. Exact results 
obtained via second order perturbation theory have been 
obtained for models with or without a cosmological con- 
stant. It is then crucial to know whether these results 
would be affected in case of quintessential dark energy. 

The aim of this paper is therefore to examine the 
growth of structure in both the linear and the nonlin- 
ear regime. To illustrate our results and their robustness 
we consider various realistic models of Quintessence. So 
far the evolution of large scale structure has only been 
studied with effective quintessence p] , ^l| , ^2[ and we will 
see that it does not provide a realistic account of what is 
happening in explicit models of quintessence. 

This paper is divided as follows, in Section II we de- 
scribe the models we use and in particular the evolution 
of the vacuum equation of state they imply. In Section 
III the results for the linear and second order growth rate 
are presented. Implications of these results are discussed 
in section IV for the nonlinear power spectrum. 



]17| and the model developed by Brax and Martin in 
which the potential shape incorporates generic Super- 
Gravity factors (l^Jl^]. In this section we succinctly re- 
view the properties of the effective, Ratra-Peebles and 
SUGRA quintessence models and compute the resulting 
equation of state of the Universe in these models. 



A. Effective Quintessence 

Models of "effective" quintessence are the simplest im- 
plementations of a non-trivial vacuum equation of state. 
It is simply assumed that the equation of state param- 
eter is fixed and represent an average value of a cosmic 
component following a complex evolution (whether it is a 
quintessence field or not). Considered as a simplified ver- 
sion of a quintessence model, this is a valid approach if, 
for some reasons, the kinetic energy and the potential are 
almost constant and of the same order. This condition 
(which is not the slow-roll condition where the kinetic en- 
ergy is much smaller than the potential) seems unlikely 
to be verified in a realistic framework. We will nonethe- 
less compare realistic models with this approximation to 
show its impact on observed quantities. 



B. Tracking Quintessence 



II. THE QUINTESSENCE MODELS 

We postulate that the content of the universe includes 
a scalar field, Q, of potential V. This scalar field is re- 
sponsible for the dark energy we observe today and usu- 
ally described as a cosmological constant. Its motion 
equation is given by the Klein-Gordon equation 



dV 

Q + 3HQ = - — 
oQ 



(1) 



and it contributes to the energy and pressure terms with 

Pq = V(Q) + \q\ 

PQ = ~V(Q) + \q 2 - 
The equation of state of the dark energy 



Pq = u Q pq 



(2) 



(3) 



is a priori no longer characterized by a constant ujq = — I 
parameter. It can vary from ujq = — 1 when the dynamic 
of the field is dominated by its potential, to loq = 1, when 
the kinetic energy dominates. In all the models we will 
consider, the parameters will be chosen so that Slo = 0.3 
and Qq — 0.7 today, except otherwise mentioned. 

In the following we focus our analysis on two models 
with tracking solutions that provide explicit time depen- 
dency of the equation of state, the Ratra-Peebles model 



In a very wide class of quintessence models the field 
dynamics exhibits a tracking solution. It is such that the 
evolution of the dark energy, during radiation and mat- 
ter domination, is completely determined by the poten- 
tial shape regardless of the initial conditions^]. In other 
word, the only tuning they require to reproduce today's 
observations is the energy scale of the potential. Eventu- 
ally, this scale will have to be explained by high energy 
physics computations. While this task seems insupera- 
ble in case of a pure cosmological constant, it mi ght be 
within theoretical grasp for tracking quintessence [|13 33 

The phenomenological properties of such class of mod- 
els can be summarized through the time evolution of the 
cosmic equation of state. The detailed behavior of the 
field in the first stages of its evolution depends on the 
initial conditions. If initially pq represents a fair fraction 
of the cosmic energy density, the field Q rolls quickly 
down its potential so that the quintessence energy den- 
sity is purely kinetic. It is slowed by the expansion until 
it freezes at a value larger than the one corresponding 
to the tracking solution. The value of the field remains 
then constant — the quintessence energy density is purely 
potential — until it catches with the attractor solution. 
Once on the attractor solution, the equation of state 
parameter of the quintessence field takes a value that 



*over hundred orders of magnitude 
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depends only on the shape of the potential and on the 
equation of state of the dominant specie of the universe 
(it therefore changes at equivalence). When the energy 
density of the field starts to dominate, the field follows 
an inflationary type slow roll solution whose equation of 
state is closing to loq = —1. These behaviors are dis- 
played on Fig. [l] for the potentials we adopted. The time 
at which the tracking solution is reached is completely ar- 
bitrary and has no effects on the quantities we consider 
in the following. 

Our analysis will be done for two tracking models: 



• the Ratra-Peebles model jlTJ whose potential is 

M 4+ a 



Vrp(Q) 



(4) 



and a = 11 which both lead to equation of state 
uiq ~ —0.8 at zero redshift. These choices of pa- 
rameter lead to an energy scale M from 10 6 to 10 11 
GeV, which does not contradict our knowledge of 
high energy physics. 

The two models have the same tracking solution and the 
equation of state parameter is thus the same, given by Eq. 
(||), on it. Differences between two models arise when the 
field leaves the tracking solution. At this time, the field 
value is of the order of the Planck mass, and the SUGRA 
correction of the latter models starts to dominate. This 
SUGRA correction cures the problems encountered by 
the Ratra-Peebles potential by quickly slowing the field 
as it rolls down thus providing a smaller equation of state 
parameter 16 ll|. 



and which is the simplest model exhibiting a track- 
ing solution. In particular, it is very hard, with this 
potential, to get an equation of state loq < —0.7 
while keeping a reasonable (from the high energy 
physics point of view) energy normalization for M 
if 17a = 0.7 today. Note that for such a poten- 
tial the vacuum equation of state of the attractor 
solution is given by 



—2 + au>B 
a + 2 



(5) 



where lob is the equation of state parameter of the 
background fluid (1/3 for a radiation dominated 
universe, for a matter dominated universe). In 
the following we will consider the case a = 2, which 
gives loq ~ —0.6 today, marginally consistent with 
the supernovae observations although it leads to an 
unrealistic low energy scale for M. 



• the SUGRA model, proposed by Ph. 
Martin Jll|[l5| whose potential is 



M 4+a 

^Sugra(Q) = -qET ex P 



4tt 



^Planck 



Brax and J. 



(6) 



The corrective factor is motivated by the fact that, 
in the Ratra-Peebles scenario, the field naturally 
reaches the Planck scale at low redshift. If the 
quintessence potential is to be derived from models 
beyond the standard model of particle physics that 
are expected to include super-gravity properties, it 
is natural to expect super-gravity corrections in the 
shape of the potential. The potential Ph. Brax 
and J. Martin proposed is actually the extension of 
the Ratra-Peebles potential, with a generic Super- 
Gravity correction (the exponential term). This 
last model is of particular interest since its pre- 
dictions are in good agreement, for a wide range 
of parameters, with the SNIa measurements. We 
studied here two examples of this potential, a = 6 
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FIG. 1. The evolution of the vacuum equation of state as a 
function of the expansion parameter a for different cosmologi- 
cal models. The dotted line corresponds to a vacuum equation 
of state, p — — 0.8p; the dot-dashed line to a Ratra-Peebles 
solution with a = 2; the dashed line to a SUGRA behavior 
with a = 6 and the solid line to a = 11. The amplitude of the 
quintessence potentials is such that ^matter = 0-3 at z = 
in all cases. 



C. Solution of the equations with Quintessence 

The quintessence field contributes to the Friedman 
equations, and therefore to the evolution of the expansion 
rate of the universe, 



" 1 3M Planck 
a An 

a 3M Planck 



Ptot. ( 7 ) 
(Ptot. + 3 Ptot.) ( 8 ) 



where p^ot is the total energy density of the Universe 
and Ptot ^ s pressure assuming we live in a zero curvature 
universe. It is convenient to define the parameter u) as 
the effective equation of state parameter of the ensemble 
of the cosmic fluids, 



Ptot. = ^Ptot. 



(9) 



3 



This parameter is expected to vary from 1/3 in the 
radiation dominated era, ui — in the matter dominated 
era tow-> —1 when the vacuum energy dominates. The 
shape of this transition, and its implication on the growth 
of structure is precisely what we investigate in this paper. 




FIG. 2. The evolution of the global cosmic equation of state 
as a function of the expansion parameter a for different cos- 
mological models for a redshift range, z = to 10. Same 
convention as Fig. [I] with thin solid line corresponding to a 
pure cosmological constant such that Qa = 0.7 at z = 0. 

The evolutions of the equation of state of the uni- 
verse are shown on Figs |l|||. They show that the tran- 
sition from a matter dominated universe to = to a 
vacuum dominated universe is much smoother in case 
of Quintessence models. In fact, the universe leaves the 
lu = line much sooner in the tracking quintessence mod- 
els than in the effective quintessence, or the A models. It 
is then natural to expect significant effects on the angu- 
lar distances or on the growth of fluctuations. Moreover, 
the low redshift behavior of the global cosmic equation 
of state is very different in the three quintessence mod- 
els. This also should induce significant phenomenological 
differences between the A models and the quintessence 
models. 

The implication of these behaviors for the angular dis- 
tances is shown on Fig. [| The differences seem not very 
noticeable at small redshift. However, they build up to 
be significant when the effect is integrated to the last 
scattering surface. The end values of the angular dis- 
tances for the different tracking quintessence scenarios 
are very close to each other, although they correspond 
to different values of the equation of state today (see fig 
||). Not surprisingly quantities sensitive to the angular 
distance at high z, such as the position of the first acous- 
tic peak, have been found to depend upon the vacuum 
equation of state p3j . One should also expect signif- 
icant differences in the amplitude of the lens effect on 
the CMB anisotropics (as its depends on the angular dis- 
tances between the observed objects and the lenses). We 
note however that a A model with a lower value of Qq 
can reproduce fairly well the low redshift behaviour of 
the angular distances in our Quintessence models. 
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FIG. 3. Evolution of the comoving angular distance in the 
different scenarios compared to a pure cosmological constant 
from 2 = 1000 to the local universe. Same conventions as 
in Figs. 1-2 are used for the line styles. The lower solid 
line corresponds to a pure cosmological constant model with 
Q.0 = 0.4. 

III. THE GROWTH OF STRUCTURE 
A. The linear growth rate 

In the previous section we have observed that the evo- 
lution of the universe with a quintessence component be- 
ing smoother than with a pure cosmological constant, its 
departure from the case of an Einstein-de Sitter (EdS) 
universe occurs later in the former. This effect has 
been described before p2| , [Tl 12 however only in the con- 
text of effective quintessence, but it is clearly amplified 
here because the energy fraction of the quintessence field 
can be much larger at high redshift in cases of realistic 
quintessence models. 

In this section we investigate the impact of these ef- 
fects upon the evolution of large scale structure. We 
only consider the large scale structure history after re- 
combination, a time at which the dark matter fluctua- 
tions dominate. After recombination and at sub-horizon 
scales the quintessence field perturbations correspond to 
decaying modes and can therefore be ignored. Within 
these assumptions the growth rate of the density con- 
trast at linear order is driven by the equation J34J, 



Z)i(t) + 2HD 1 {t) - -H 2 Q(t) Dx(t) = 



(10) 



where il(t) corresponds to the fraction of energy in the 
matter component. The growth rate is independent on 
the wavelength of the fluctuations as long as we con- 
sider fluctuation at sub-horizon scale and if we neglect 
the pressure effects. The time evolution of Di provides 
the amplitude of the density fluctuations. Fig. || gives 
the growth rate, for different models, as a function of 
redshift compared to the growth rate in EdS[| In mod- 



^ whose solution is well known -Di(t) = a(t) 
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computation of higher order terms in a perturbative ap- 
proach. More precisely the reduced skewness defined as, 
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FIG. 4. The ratio D + /a (normalized to unity at z > 100) 
for the different scenarios. The today growth rate is smaller 
by about 20% in tracking quintessence scenarios. 



els with a pure cosmological constant, the growth fac- 
tor remains very close to the EdS solution for a long 
period, then changes abruptly between redshift 2 and 
3. In quintessence models the solutions follow the same 
scheme, yet with a smoother slope change. However, for 
the same normalization at z > 100, when the Universe 
is very close to a Einstein-de Sitter model, the today's 
growth rates are quite different. The quintessence sce- 
narios with a tracking field exhibit a smaller growth to- 
day (of order 20 to 30 percent less). To say it in other 
words, for the same trg, the quintessence models demand 
for larger density fluctuations at early times. This effect 
should lead to a difference between CMB normalizations 
and low redshift normalizations. 

The origin of this difference is clear. It is due to the fact 
that the energy fraction in the quintessence field remains 
significant for a much longer time. With this respect the 
effective quintessence solution reveals very similar to the 
A scenario, whereas realistic models of quintessence leads 
to linear growth rates that depart from the EdS case 
at redshift as large as 30! Clearly, models of effective 
quintessence that can match the SNIa observations do 
not provide a good account of the linear growth rates 
found in realistic models of quintessence. 



B. Second order growth rate 

As mentioned in the introduction, for Gaussian ini- 
tial conditions the second order growth rate determines 
the rate at which non-Gaussian properties emerge in the 
matter density field. In a perturbation theory approach 
it is indeed this quantity that determines the value of 
the large-scale skewness. Furthermore, it turns out that 
weak lensing surveys can be used as a test-ground for 
this effect and would provide a robust constraint on Qq 
through the value of convergence skewness We 
examine here to what extent this approach remains valid 
in a quintessence cosmology. 

Such perturbation theory calculations are based on the 



(11) 



can be related to the second order growth rate |36| and 
more specifically to the second order growth rate in the 
spherical collapse dynamics |3?],[38| . In this case it is sim- 
ple to expand the local density contrast to second order 
with respect to the initial density fluctuations, 



6sc(t) = D 1 (t)6 l + ^6* + 



(12) 



with time dependent coefficient that can be explicitly cal- 
culated for any cosmological model. The function D 2 (t) 
is the growing mode of the equation, 

D 2 (t) + 2HD 2 (t) - ^H 2 n(i) D 2 (t) = 

2,H 2 n(t)Dl{t) + \bl{t). (13) 

The 3D density skewness at large scale (and when 
smoothing effects are neglected) is then directly propor- 
tional to D 2 (t), and is given by 



S* =3 



D2(t) 

my 



(14) 
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FIG. 5. The skewness for the different scenarios. The hor- 
izontal straight line gives the skewness in EdS models. The 
variations on the quintessence models are very small. 

From Fig. || it is clear that the variations of the skew- 
ness with the cosmological models are very small (below 
percent level) and are likely to remain undetectable. It 
means that the second order growth rate does not in- 
troduce further dependence with the vacuum equation of 
state. This result actually extends a property already 
known for the dependence of 5*3 on Slo for open universes 
|p6fl , flat universes with a cosmological constant (3|] or 
in some flavors of non standard vacuum equation of state 
& 

The skewness of the convergence is thus expected to 
be left unchanged, except through the dependence of the 
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angular distances and the linear growth rate of the fluc- 
tuations. We recall here the formal expression of the 
convergence skewness in perturbation theory for a power 
law spectrum (of index n) pfl]. For a flat Universe, it is 
given by, 



S3 



f" s dVw 3 (V) D\{V) 2?- 2 («+ 2 ) [Sj D - |(n + 2) 



tf-dVv?(V)Dl(V)V- 



(n+2) 



with an efficiency function w(T>) given by 

V) 



w(V) 



3 n V{V S 



(15) 



(16) 



where T> s is the comoving distance to the sources and 
S3 is the skewness parameter for the 2D dynamics. The 
latter can be related to the 3D one, since it implies only a 
different combination of the terms appearing for the 3D 
case p8|,p0i. It finally gives, 
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s 3 = 77 + 7 



4Df' 



(17) 



which amounts to 36/7 for an Einstein-de Sitter case. 

TABLE I. Value of the skewness, Eq. (|l|) of the local 
convergence in weak lensing surveys for sources at redshifts 1 
or 2 and for a power law index n = —1.5. 



skewness 



Z s — 1 



z s = 2 



A model (fi = 0.3) 
A model (fi = 0.25) 

w Q = -0.8 
Ratra-Peebles, a — 2 
Sugra, a = 6 
Sugra, a = 11 



76. 

85. 
83. 
91. 
85. 
86. 



26. 
28. 
28. 
32. 
30. 
30. 



In table | we present the expected skewness for the dif- 
ferent models we have considered for sources at redshifts 
1 or 2. The results show that the projection effects on the 
value of the skewness can be pretty large. They increase 
the value of the skewness so that Quintessence models 
with Qq = 0.3 mimic what one expects for a fio = 0.25 
model with a pure cosmological constant. 

It is to be noted that angular diameter distances, in a 
quintessence scenario, rather resemble a A model with 
a larger value of Qq (see Fig. |||). From those two 
joint observations it should then be possible to test the 
quintessence model hypothesis. Results should however 
be extended to the intermediate and nonlinear regime 
where most of the data are going to be although we ex- 
pect the qualitative results found here to remain valid. 



IV. NON-LINEAR MATTER POWER 
SPECTRUM IN QUINTESSENCE MODELS 



A. The shape of the nonlinear power spectrum 

We complete these investigations with the non-linear 
evolution of the matter density contrast power spectrum. 
Assuming stable clustering ansatz and the Hamilton et 
al mapping ETJ], we compare the power spectrum in 
quintessence and cosmological constant models in the 
deeply non-linear regime. 

The linear matter power spectrum in quintessence sce- 
narios had been studied before {^@||j23). For modes 
inside the horizons, models with a pure cosmological con- 
stant or with a quintessence field show very little differ- 
ences in the shape of the linear transfer function. Lim- 
iting our study to those modes, we can reliably approxi- 
mate the linear quintessence power spectrum by a stan- 
dard ACDM one. 

However there is no reason for the non-linear evolu- 
tion of these models to lead to the sa me power spectrum. 
Indeed, we showed in section [II A that generically the 
large scale structure grows more slowly in quintessence 
scenarios. Hence, for the same amount of structure to- 
day, the density contrast had to be bigger in quintessence 
scenario at early time. This implies that modes that are 
in the non-linear regime now have reached this regime 
sooner in quintessence scenarios. 

Following the idea of Hamilton et al. j|0|, later ex- 
tended by Peacock and Dodds (4^j44|], we postulate 
that one can describe the effects of non-linear evolution 
through a universal function / n i that maps the linear 
power spectrum onto the non-linear one 



-1/3 



A nl (fc n l) = /nl(A 2 (fc)) 

fc=(l + A nl (fc nl )j fc nl 
A 2 (fc) = 47rfc 3 P(fc). 



(18) 



Enforcing stable clustering, the previous authors have 
shown that this function must follow an asymptotic be- 
havior at large x such that f n \(x) oc g($V)~ 3 x 3 / 2 — where 
g(£l) — Di(a)/a, is the ratio of the linear grow th fac tor 
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to the EdS growth factor described in section 
and have proposed analytic forms for / j that depend 
on the cosmological parameters through g(a) only and 
that are calibrated on various A^-body simulations. We 
assume that their results hold for quintessence scenarios. 
In particular, we assume that the non-linear regime al- 
ways reaches a stable clustering regime. Moreover, and 
in absence of quintessence N-body simulation for our par- 
ticular scenarios, we also assume that the normalization 
factor in the asymptotic branch is independent on the 
cosmological scenario. These assumptions are not triv- 
ial and can probably be challenged (see for instance J32j l 
where the behavior of the nonlinear power spectrum is 
investigated for effective quintessence models with a dif- 
ferent perspective). 

From the Hamilton et al. ansatz we expect very dif- 
ferent behaviors for the small scale non-linear power 



6 



spectrum when quintessence and non-quintessence mod- 
els are compared. We previously obtained that about 
20% to 30% discrepancy is expected, depending of the 
quintessence potential, between g® and D A . The conse- 
quences of this a priori modest discrepancy are dramatic 
for the nonlinear power spectrum. For a mode that en- 
tered the non-linear region long before redshift z, we have 



A m (fcnl>*) 2 
gQ(z)-z AQ(k Q ,zf 



(19) 



S A( z )-3 A A (k A ,zf 

where /cq (resp. k\) is the linear mode in the 
quintessence (resp. A) linear power spectrum giving rise 
to the fc m mode in the non-linear spectrum. Assuming, 
as stated before, that at the sub-horizon scales we are 
interested in, the linear power spectrum of the models 
are identical, up to a normalization factor Pq, we write 
P(k, z) = g 2 {z)a 2 Pq k n with n > —3. For a mode in the 
non-linear region we have 



a 2 P 



(20) 



and Eq. (19) reduces to 



Q( Z y 3 g^(z) 2 a 2 P^k^ 



Qj.n+3\3/2 



^nl^nl' 



9 A (z)- 3 \g A (z) 2 a 2 P A kl +3 
. S.( 1 _n+3\ 

Q\ 2V 1 n +5) 




(21) 



We suppose here that the spectral index n is identical 
for both P < 3(/sq) and P a (/ca) which is a reasonable ap- 
proximation. Finally if we set the spectrums to both fit 
the cluster normalization, the ratio Pg / P^ is simply the 
ratio of the growing modes at z — and we get 



^>nb 



P^k^z)' \g A (z = 0) 



g Q (z = Q) 



-3(i-5S) 



(22) 



Note that the exponent gets close to —3 when the spectral 
index goes to —3. Given the variation of the ratio 
given in Eq. ( p2|) can be as large as 2! 

The same ratio is easier to compute for modes in the 
linear regime, 



pQ(k,z) fg Q {z)g A (0) 



P A (k,z) \gQ(0)g A (z) 



(23) 



These simple investigations show unambiguously and 
with a limited number of assumptions that the shape 
of the nonlinear power spectrum is very sensitive to the 
presence of a quintessence field. In order to have a full 



description of the power spectrum behavior, including 
the intermediate regime, we use the Peacock and Dodds 
prescription. This formula has been shown to be reason- 
ably accurate for effective quintessence models in N-body 
simulation, at least at low redshift [p2| . Anyway, we are 
more interested, in this paper, in the general trend rather 
than a percent accurate description of the transition be- 
tween the linear and non-linear evolution. Fig.^j shows 
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FIG. 6. the ratio P® D /P£ D at * = 0. The solid and dashed 
lines are the SUGRA quintessence models for a = 11 and 
a = 6, the dot-dashed one is Ratra-Peebles and the thin dot- 
ted line is the effective quintessence model for ujq = —0.8. 
We also show here the results for a Qo = 0.4 flat ACDM 
(thin dashed line) and the O = 0.25 flat ACDM (thin dou- 
ble dashed line) models which where discussed in previous 
sections. 

a comparison between quintessence models and a stan- 
dard ACDM. The curves represent the ratios between a 
quintessence non-linear P 3 £> at z — and the ACDM 
one assuming the linear power spectra have the same 
normalization. As expected, the small scale ratio tends 
toward the (g+/g A )~ 3 asymptote. The shape of the curve 
in this region is described by the 3 (1 — (n + 3)/(n + 5)) 
power calculated above. Since the effect goes like the 
growth ratio to the third, it is much smaller for effective 
quintessence where the discrepancy between its growth 
ratio and the A model growth are smaller. We stress 
therefore again that the use of a constant equation of 
state cannot account of the amplitude of this effect as it 
is expected in realistic models. 

Note that the transition between the linear and non- 
linear regime depends on the ansatz used for / m . How- 
ever, a very sharp transition, as encountered here, is not 
unnatural. It accounts for different times a given mode 
enters the non-linear regime in different models. If one 
follows a given mode throughout its evolution, it will first 
obey the linear growth and evolves as a 2 g 2 (a). Then it 
enters the nonlinear regime and grows as a 3 . The tran- 
sition between this two regimes is very sharp. When 
models with a different growth factor are compared, this 
rapid transition translates into a sharp increase of the 
power spectrum ratio between the linear and the non- 
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linear regime. 

Moreover we note that this effect cannot be mistaken 
with a variation of Hq. The latter has a much more dra- 
matic effect on the shape of the linear power spectrum 
through a change of the shape of the transfer function. In 
this case, not only is the linear growth rate changed but 
the position of the maximum of the linear power spec- 
trum is also shifted. 

In principle large-scale galaxy surveys such as the 2dF 
or the SDSS should be able to put constraints on the am- 
plitude and shape of the power spectrum. The possible 
effects of biasing mechanisms, that are extremely badly 
understood in the transition regime between the linear 
and the non-linear regime, prevent however to have a 
robust and reliable test of these scenarios. In the next 
section we rather try to validate these properties in the 
context of a better defined observational procedure: the 
projected power spectrum in weak lensing surveys. 



B. Projected power spectrum in weak lensing 
surveys 

Weak lensing surveys can potentially provide us with 
precision maps of the projected density up to redshifts 
around one [1451 [MEG 49 27Lp(J] . These measurements are 



expected to be free of observation biases once the redshift 
distribution of the sources is known. 

Weak lensing surveys, through the observation of the 
deformation of background galaxies, can give access to 
the convergence field. The latter can be written as the 
projection of the matter density along the line-of-sight 



k [a) 



dVw(V)6(a,z(V)), 



(24) 



where the V stands for the angular distances (and this 
formula is valid for a flat spatial curvature only). The 
geometrical kernel w(T>), defined in eq. (46), accounts 
for the projection effects. 

We expect that the ratio of the power spectrum of the 
convergence field for models with a cosmological constant 
and models with a quintessence field exhibit roughly the 
same properties than the ratio of the three dimensional 
power spectra. However, the value of this ratio is ex- 
pected to be affected by a corrective factor induced by 
the geometrical kernel that itself depend on the details 
of the cosmological model. In the following we present 
an evaluation of this rescaling factor along the line of 
reasoning of the previous subsection. 

In the small angles limit, we have 



dVw(V) 2 lAl_P 3D ( k ,z(V))e £vs (25) 



h J ( 27r ) 

so that the convergence power spectrum is simply 



dX> 

V 2 



w(V) 2 P 3D (l/V,z{V)), 



(26) 



where a possible redshift evolution of the shape of the 
power spectrum is included. The kernel w is a bell shaped 
window that reaches its maximum at z e g = z(T> s /2). To 
evaluate roughly the rescaling factor we will approximate 



?/>• 



by a simple Dirac function w (T>) 



w 



eff 



5{V-V s /2) 



with w e g = J 3 dT>w{T>) and z e g- ~ 0.4 (it depends on 
the cosmology we are considering) for sources at redshift 
z s . Now, the ratio of the convergence power spectrum 
simply writes, 



w 



Q 

Off' 



/Vf\ P^ D {21 



z cff 



"eff 



Pf D (2l/V$,z A 



eff 



(27) 



For a mode in the non-linear region, a few percent error 
in the position of the mode is not significant, so that we 
can ignore the difference between Dq and D\ in the last 
term of the equation. As a result 




9 Q (z = 0) 

g A (z = o) 



-3(1-58) 



(28) 



which, compared to Eq. (£2|), contains an extra geometri- 
cal factor due to the projection effects. It evaluates to 0.8 
to 0.9 depending on the models and position of the source 
plane. We expect therefore the conclusions reached in the 
previous section to survive in weak lensing observations. 
Similarly, for a mode in the linear region, we have 



i?(0 



w 



off 



„ 5 A (0) 



(29) 



^ gQ(Q) . 9 A(zA ff) J ■ 

Table || gives the value of the expected ratio in the linear 
region. It is about 0.9. It indicates that the normaliza- 
tion ratio for the linear 3D power spectrum and the one 
for the projected weak lensing spectrum will differ by this 
amount. 

TABLE II. Evaluation of the ratio P$ / 'P A in the linear 
domain, from Eq. (H). 





2 = 1 


z = 2 


z = 1000 


ujq = -0.8 


0.94 


0.95 


1.06 


Ratra-Peebles a = 2 


0.86 


0.88 


1.20 


SUGRA a = 6 


0.92 


0.92 


1.16 


SUGRA a = 11 


0.91 


0.91 


1.18 



These semi-analytical results give a good account of 
the ratio P^/P^ for all angular scales. At large scale, 
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it is roughly flat; its value is given in table H Then, 
as we get closer to the transition between the linear and 
the non-linear regime, the ratio exhibits a slight drop. 
Indeed, from Eq. (^7j) the ACDM model enters the non- 
linear regime before, because of the difference between 
the z e g. Hence, one expects to have, for a few modes, a 
P£ that rises quicker than its quintessence counterpart. 
Then, when the quintessence power spectrum also hits 
the non-linear regime, the ratio will exhibit a shape very 
similar to the one of the three dimensional power spec- 
trum ratio.with the rescaling factor computed above. 

In Fig. [7| we present the explicit computation of the 
nonlinear power spectra of the convergence field using the 
prescription of Peacock & Dodds to compute the redshift 
evolution of the 3D power spectrum. Unlike Fig. ^, the 
power spectra are not cluster normalized. In this case 
the power spectra are normalized so that weak lensing 
amplitudes match at 10' scale when computed with the 
linear power spectrum and match the amplitude of the 
recent detections of weak lensing effects (e.g. as ~ 1 for a 
A-CDM model with A = 0.7). Because of the projection 
effects given in Table 0, this is not equivalent to normal- 
ized linear 3D power spectrum. Projection effects also 
slightly change the shape of the projected linear power 
spectrum. The redshift of the sources is simply here as- 
sumed to be unity. The differences in the shape of the 
power spectra is clearly visible and should be already 
within observational constraints. 

We also give, following the same prescription for the 
normalization, the effects of a change of Oo- In this case, 
because we normalized to the convergence linear power 
spectrum, the effects of Slo also directly affects the nor- 
malization. Compared to 3D power spectra, it actually 
worsen the situation and make the distinction between 
quintessence models and such models striking. 




The result of fig. gives us hope to constraint strongly 
the quintessence scenario using weak lensing surveys. 
The next generation weak lensing surveys will made avail- 
able wide surveys where a precise determination of the 
lensing effect will be possible for a range of scale large 
enough to map the sharp rise predicted here. For exam- 
ple, measurement of the weak lensing effect amplitude 
at one degree scale and at one minute scale with only a 
ten percent precision appear sufficient to test a SUGRA 
quintessence hypothesis. It seems that the observational 
requirements are much more modest than direct mea- 
surements of the angular distances through SNIa obser- 
vations. 



C. Weak lensing on Cosmic Microwave Background 

In passing we note that weak lensing effects on CMB 
maps could be used also to test quintessence hypothesis. 
Amplitude of the effects are mainly given by the ampli- 
tude of the fluctuations of k, a 2 , along the line-of-sight 
pllp2|. In Table III we show the amplitude of the lens 



effects on the last scattering surface at two different an- 
gular resolutions. They are mainly sensitive to the linear 
change of the growth rate integrated over the line of sight. 
It would probably be not a crucial test for the nature of 
the vacuum energy but it is potentially an important test 
to pass once the general cosmological parameters will be 
determined. If the coming generation of observations call 
for quintessence, observation of an excess of power of the 
lens effect as suggested by these calculations, would be 
an important consistency test. 

TABLE III. Ratios between the amplitudes of the lens ef- 
fects on the last scattering surface for different models and 
the standard ACDM model (A = 0.7) at two different angu- 
lar resolutions. 



at z 



1000 



10' 



fi = 0.25 A = 0.75 model 
n = 0.4 A = 0.6 model 
uq = -0.8 
Ratra-Peebles a — 2 
SUGRA a = 6 
SUGRA a = 11 



1.23 
0.68 
1.20 
1.49 
1.29 
1.33 



1.29 
0.73 
1.21 
1.54 
1.32 
1.36 



1000 10000 100000 



V. CONCLUSION 



FIG. 7. Ratios P$(l)/P£(l) for a source plane at z = 1. All 
models are normalized so that a K are the same at 10' scale 
and correspond to a Slo = 0.3, flat ACDM with as = 1. The 
observational window, in which measurements with an accu- 
racy better than 10% is foreseeable, corresponds to £ of 200 
to 10,000 (minute to degree scale) where the most dramatic 
changes take place. 



In this paper we have examined the growth of structure 
in quintessence models in both the linear and the second 
order regime and present their more striking implications 
for the statistical properties of the low redshift large-scale 
structure of the universe. 

We paid particular attention to cases of realistic imple- 
mentations of quintessence field since they lead to scenar- 
ios where the energy fraction in the quintessence compo- 







nent can represent a significant fraction of the total en- 
ergy density over a long period. We indeed found that 
this effect is responsible of important differences in the 
behavior of the linear growth rate of the fluctuations: 
For the same values of fio, realistic quintessence models 
lead to a linear growth rate that can be 20 or 30% lower 
compared to models with a pure cosmological constant 
or with an effective quintessence component (where the 
vacuum has a constant equation of state which matches 
the angular distances constraints). 

Consequences of this discrepancy have been explored 
at the level of the nonlinear power spectrum for which 
such differences are amplified. For power spectra with 
identical linear normalization (at z — 0), the variation of 
the amplitude of the nonlinear power spectrum can be as 
large as 2. 

We have also computed the second order growth rate 
of the fluctuation. We found that, when expressed in 
terms of the square of the linear rate, it is not sensitive 
to the nature of the dark energy. This ratio is actually 
not significantly sensitive to any of the cosmological pa- 
rameters. In this respect our result extends previously 
known properties. 

Weak lensing surveys appear to be the natural play- 
ground for such effects. They combine effects on the an- 
gular distances and on the growth rate of the fluctuations. 
We show that the skewness of the convergence field, at 
large angular scale, is notably sensitive to the projection 
effects. It is to be noted however that a universe with 
quintessence field does not resemble a universe with a 
cosmological constant and larger matter density (as it is 
the case for the behavior of the angular distances) but 
rather with a lower density parameter. 

Moreover the shape of the power spectrum of the con- 
vergence field, which identifies with a projected 3D mat- 
ter power spectrum, retains the properties found for the 
3D nonlinear spectrum. It appears clearly that CDM 
family models for a flat universe can be distinguished 
from one another: a variation of f^o changes the shape 
of the linear power spectrum, whereas the introduction 
of a quintessence field changes the time at which modes 
become nonlinear. 

The precision level of the current semi-analytical pre- 
dictions for the shape of the nonlinear spectrum does not 
permit so far to make precise predictions from which the 
quintessence potential could be reconstructed. Moreover 
the use of the prescription of Peacock & Dodds for models 
of quintessence with a tracking solution should probably 
be validated with specific numerical simulations. 

The calculations have been done in two specific models 
of quintessence, the Ratra-Peebles model and the Sugra 
model developed in [fl6| . We think however that our con- 
clusions would survive for any model where the energy 
density in the quintessence component can be a signifi- 
cant fraction of the total energy up to recombination. 
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